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۲ (نویسنده مسئول) 


حکیده 


در مطالعه حاضر به<"وسیله یک مدل ت رکیبی نحوه انتقال رسوبات درون مجاری فاضلاب "رو مد ل-سازی م ی-شود. به‌بیان دیکس 
سه‌فازی (هواه آب و رسوب) ارائه داده می"شود. در این مطالعه, از الگوریتم ژنتیک برای افزایش کارایی سیستتم استنباط فازی 
عصبی تطبیفی با تنظیم نوابع عضویت و سپس به حداقل رساندن مقادیر حطا استفاده گردیده است. الکوریتم تیک (0۸) یک 
روش فرابتکاری با الهام از فرایند انتخاب طبیعی می-باشد که متعلق به طبقه-بند ی گستردهستری از الگوریتم-های تکاملی (5۸) 
است. سپس با استفاده از مقادیر ورودی» ۱۲۷ مدل ت رکیبی ۸۷5-0۸ معرف یگردید. همچنین مقادیر درصد میانکین مطلق عطا و 


تعطای جذر میانگین مربعات برای مدل برتر بهترتیب مساوی ۵/۵۲۹ و ۰/۳۱۵ بدست آمد. 


واژه [ آهای کلیدی: رسوب. کانال دایروی. انفیس. الگوریتم ژنتیک» مدل ترکیبی 
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معد مه 

به‌منظور استفاده از معیار انتقال رسوب در حالت حد ته- 
نشینی. نیاز به اطلاعاتی مانند اندازه دانه‌هاه چگالی» غلظت 
رسوب و عمق جریان مورد نیاز است که بتوان یک معادله 
برای پیش‌بینی نرخ انتقال رسوب به‌عنوان بار بستر یا بار 
معلق ارائه داد. انتقال رسوب در حالت حد ته‌نشینی به دو 
صورت بار معلق و بار بستر می‌باشد. به‌منظور استفاده از 
معادلات خودشویی موجود. حالت انتقال (بار بستر با بار 
معلق) باید مشخص باشد. کلی و گولارته یک مطالعه 
آزمایشگاهی در ارتباط با فرسایش بستر کانال‌های رسوبی 
اجرا نمودند. آنها بیان نمودند که با تغییرات عکس دما 
مقدار فرسایش بستر رسوبی کاهش می‌یابد[ ۱ راودکیف و 
تان در یک تحقیق آزمایشگاهی نحوه فرسایش رسوبات 
چسبنده درون کانال‌های دایروی را مورد ارزیابی قرار دادند. 
آنها تغییرات فرسایش رسوبات را در مقابل ۳۳ جریان برای 
شرایط مختلف اندازه‌گیری نمودند[ ۲]. کرابتری در یک 
مطالعه آزمایشگاهی کانال‌های فاضلاب را به پنج نوع 
مختلف بر حسب بستر تقسیم‌بندی نمود[۳]. آلوارز 
آزمایش‌های انتقال با رسوبات غیرچسبنده روی بسترهای 
زبر و صاف در لوله‌ای به قطر ۱۵۴ میلیمتر با جریان قسمتی 
پر انجام داد. او برای داده‌های بستر صاف. رابطه‌ای که 
بهترین برازش با کاربرد آنالیز رگرسیون حطی چندگانه 
داشته باشد را با استفاده از آنالیز ابعادی بدست آورد[۴]. 
عظمت ا. و همکاران [۵] با استفاده از ترکیب داده‌های آب 
قانی [۶] و وونگ ویسس امجای و همکاران [۷] و استفاده 
از آنالیز رگرسیون, رابطه ارائه شده توسط آب قانی را 
اصلاح کردند. اتا و پروسکیا تتش برشی بی‌بعد به‌منظور در 
نظر گرفتن رسوبات با اندازه 4 که روی بستر زبر با زبری 1 
(0>) حرکت می‌کنند. آنها در مطالعه خود تغییرات دبی 
جریان را بر روی نرخ انتقال رسوبات را اندازه‌گیری 
تمردنلل قزرجمه اتتیر تکیک‌های گوناکزن شنگته متسین 


و هوش مصنوعی در پیش‌بینی و الگوشناسی پدیده‌های 
پیچیده مشکل و غیرخطی علوم مختلف توسط بسیاری از 
ره اه فان کی ماع ]مرن 
های نرو- فازی و الگوریتم‌های مختلف شبکه عصبی 
مصنوعی برای حل مسائل مختلف علم هیدرولوژی و 
هیدرولیک بکار گرفته شده‌اند. محققین مختلفی همچون 
شممی‌الدین (۱۹۹۷) ها کیرتولیسن فد لاوسلی (6۳۰۳۵) 
[۱۰] با استفاده از مدل شبکه عصبی مصنوعی به پیش‌بینی 
فلالهای رواتات فاز بر داعست: ات هن هی کارا 
(1۹۹۹) [۱۱] و بی و همکاران (۱۲[)۲۰۰۷] به کمک 
الگوریتم‌های مختلف شبکه عصبی جریان به داخل مخزن را 
تخمین زدند. بار رسوبات معلق روزانه با استفاده از مدل‌های 
نرو- فازی توسط کیسی (۲۰۰۵) پیش‌بینی گردید [۱۳]. در 
این مطالعه. حداقل سرعت مورد نیاز (عدد فرود) جهعت 
جلوگیری از ته‌نشینی رسوب در بستر مجاری فاضلاب‌رو به 
وسیله روش "های مختلف محاسبات نرم شامل سیستم‌های 
استنتاج فازی عصبی تطبیقی و یک روش ترکیسی مبتنی بر 
ترکیب دو روش سیستم‌های استنتاج فازی عصبی تطبیقی 
(انفیس) و الگوریتم ژنتیک برای اولین بار پیش"بینی 
می"شود. همچنین. نتایج محاسبات نرم در پیش‌بینی عدد 
فرود جریان سه‌فازی مورد نیاز جهت جلوگیری از ته‌نشینی 
رسوبات ورودی به مجاری فاضلاب‌رو با اندازه‌گیری"های 
آزمایشگاهی مقایسه می‌گردد. با مرور مطالعات گذشته‌گان 
در زمینه انتقال رسوب در مجاری فاضلاب‌رو مشسخص 
می "شود که تاکنون مدلسازی مکانیزم انتقال رسوب 
جریان سه فازی درون کانال‌های دایره‌ی شکل توسط این 
روش ترکیبی انفیس- الگوریتم ژنتیک انجام نشده است که 
در این مقاله. حداقل سرعت مورد نیاز (عدد فرود) جهعت 
جلوگیری از ته‌نشینی مواد جامد ورودی به مجاری 


فاضلاب‌رو توسط این روش ترکیبی هوش مصنوعی مدل- 


سازی می‌شود که همین بحث نشان دهنده نو اور موضوع 


مواد و روش‌ها 

ساختار سیستم استنتاج فازی عصبی تطبیقی 

در این مطالعه. به‌منظور پیش‌بینی عدد فرود مورد نیاز جهت 
جلوگیری از ته‌نشینی رسوب در بستر کانال‌های فاضلاب. از 
یک سیستم استنتاج فازی عصبی تطبیقی استفاده می‌شود. 
این سیستم به صورت یک چارچوب مدل‌سازی است که به 
صورت ترکیبی منطق فازی و شبکه عصبی مصنوعی ارائه 
شنه استه. این تال ثر کی ابه سور غلبیه کر مها 
موجود در هر دو روش منطق فازی و شبکه عصبی 
مصنوعی, ارائه شده است. در واقع با الهام از سیستم‌های 
فازی دانش اولیه در یک مجموع‌ای از قبود به منظور کاهش 
فضای جستجوی بهینه‌سازی نشان داده می‌شود. در حالی‌ که 
کیک یا فان تافقه با استقاده ای عس اتشارن آوشت که صرح 
مصنوعی الهام گرفته شده است.در این مدل از شبکه عصبی 
مصنوعی به‌منظور تنطیم کردن توابع عضویت. استفاده می- 
شود (سینگ و همکاران. ۲۰۱۲). همچنین وجود توابع 
عضویت فیرخطی در روش نرو- فازی. منجر به کاهش 
قابل توجه هزینه پیاده‌سازی یک طرح ساده بر پایه قوانین و 
حافظه مورد مصرف می‌شود. بنابراین واضح است که ترکیب 
شبکه عصبی و سیستم‌های فازی محدودیت‌های هر یک از 


این دو روش را کاهش می‌دهد و منجر به پیشنهاد یک روش 
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تعداد گره‌های لایه اول با استفاده از تعداد ورودی‌ها و تعداد 


توابع عضویت (0) برای هر یک از ورودی‌ها تعیین می‌شود. 
اما تعداد گره در ساير لایه‌ها (لایه‌های ۲ الی ۴) به تعداد 


قوانین (۴) در مجموعه قوانین فازی تیشتت رم دارد. ساختار 


داده‌کاوی برای حل مسائل پیچیده مهندسی می‌شود. سیستم 
استنتاج فازی عصبی تطبیقی» یکی از روش‌های شناخته شده 
در ترکیب همزمان شبکه عصبی و سیستم‌های فازی است. 
این روش جهت شناسایی رفتار سیستم‌های غیرخطی که با 
استفاده از مجموعه داده‌های ورودی و خروجی تعریف شده 
برای مدل مورد استفاده قرار می‌ گیرد. سیستم استنتاج فازی 
عصبی تطبیقی مدلی ساختار یافته از سیستم استتاج فازی 
است. دو تا از معروف‌ترین سیستم‌های استتاج فازی که در 
ساختار سیستم استنتاج فازی عصبی تطبیقی مورد استفاه 
قرار می‌گیرد شامل ممدانی هستند. به‌دلیل اینکه سیستم 
استنتاج فازی تاکاگی- سوگنو- کانگ نیاز به قوانین کمتری 
دارد و همچنین از روش‌های آموزش مانند پس‌انتشار که 
فکمه شاه هسسته ابسناده سکن این میت از سا کین 
بیشتری برخوردار است. به‌منظور آموزش سیستم استنتاج 
فازی عصبی تطبیقی دو روش پس‌انتشار و هیبرید که به 
صورت ترکیبی از دو روش پس‌انتشار و حداقل مربعات 
ازانه شبده است» بعنوان روش‌هاق کلاسیک نامه می 2 
شوند هر چند که امروزه استفاده از الگوریتم تکاملی به 
عنوان الگوریتم فرا ابتکاری قدرتمند در بهینه‌سازی سیستم- 
های غیرخطی. افزایش يافته است. مجموعه قوانین در نظر 
گرفته شده برای یک سیستم استتاج فازی از نوع تاکاگی- 
سوگنو- کانگ به‌صورت دو قانون اگر- آن‌گاه به‌صورت زیر 
بیان می‌شود: 

60 ,ظ ور ( 20 . ریش که ۲ 7 :1 ۴۷۱۵ 
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لایه‌های یک سیستم استتاج فازی عصبی تطبیقی به‌صورت 
زیر می‌باشد: 
لایه اول: ورودی :26 که دارای برچسب نزن می‌باشد درجه 
عضویت هر مجموعه فازی را نشان می‌دهد. تابع نود در این 


لایه را می‌توان به‌صورت زیر بیان کرد: 


,1۳۳0۷45 0 ,شرع 3 ( )ره < )0 
12,۰۰۷ < [ 
که زنل برابر با [ امین تابع عضویت برای ورودی :6 است و 
او روج گره [فاست, 
به‌دلیل عملکرد خوب و رضایت بخش تابع عضویت گوسی 
در کاربردهای مختلف مهندسی. تابع عضویت مورد استفاده 
در این مطالعه نوع گوسی انتخاب شده است. رابطه ریاضی 
این تابع عضویت به‌صورت زیر می‌باشد: 
)۳( (/(»- ۲) -)«ده < ( 2 ۸۸ 
که 2 و 0 مجموعه پارامترها هستند. 
لایه دوم: هر گره در اين لایه () که به صورت گره‌های 
دایره‌ای شکل (7) ارائه شده‌اند با استفاده از ورودی‌های 
دریافتی» خروجی را تولید می‌کند. 
( ۳ )پیک..( جک )یل( ر )رل > ,21۷ 06 
2,۱ < 6]..2 و..,2,< 1 
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لایه سوم: در این لایه گره ۳ نسبت قدرت آتش قانون *1 
به قدرت آتش کل قوانین را تعیین می‌کند. 

(۶) ,6,2 ۲۷ +..+ ۲۷+ ,۷) 0-۷-۱۷ 
لایه چهارم: هر گره در این لایه به‌عنوان یک خروجی وزن- 
دار شده سیستم استتاج فازی به‌صورت زیر عمل می‌کند: 
۳0( 2-۷ 0 
که ,۷ خروجی لایه غیرفازی‌ساز و 1 خروجی قانون 6 ام 
سیستم استناج فازی تاکاگی -سوگنو- کانگ است. برای 
تعداد 0 ام سیستم استناج فازی تاکاگی- سوگنو- کانگ؛ 
قوانین به‌صورت زیر بیان می‌شوند: 

ب( , ی۸ فف ر۱6) 0صهرل رش فا ,ک) گ 


016: (۸ 


۳۵:۸ ( بیش هف )۵00 
]<1 
12 و 1 مجموعه پارامترها هستند. پارامترهای این لایه 


تحت عنوان پارامترهای بخش تالی. شاخته می‌شوند. 


لابه پنجم: در این لایه گره دایره‌ای تک » 2 همه خروجی 
ها را به‌عنوان مجموع همه ورودی‌ها به‌صورت زیر محاسبه 
تن کت 


0 ۲ ار زا در 0 <رع‎ (٩ 
1<1 

که ۲ نشان دهنده تمام خروجی‌ها شبکه اشتتا: حهت 

شاخص میانگین مربعات خطا که به‌صورت زیر محاسبه می- 


شود استفاده می‌شود. 


(0۳ 7 ِ سس <- ۷/8 


1<ز 
که ظ تعداد نمونه‌هاه ۳۱۳۸۰ مقدار عدد فرود مشاهده شده 


در تست‌های آزمایشگاهی و 1۳۵۵4 عدد فرود پیش‌بینی 


شده با استفاده از سیستم استنتاج فاز عصبی تطبیقی است. 


الگوریتم ژنتیک 

الگوریتم ژنتیک روش قدرتمندی را برای توسعه اکتشافی 
مسائل بهینه‌سازی ترکیبی مقیاس بزرگ فراهم آورده است. 
یک الگوریتم ژنتیک مسئله را به‌صورت مجموعه ای از 
رشته‌ها که شاما رات ریرهستد کد گذاری یکت سین 
برای تحریک فرایند تکامل تدریجی. تغییراتی را روی رشته- 
ها اعمال می‌کند. در مقایسه با الگوریتم‌های جستجوی 
محلی. در جستجوی عمومی که تنها یک راه حل قابل قبول 
وجود دارد. الگوریتم‌های ژنتیک جامعه‌ای از افراد را در نظر 
می‌گیرند. کار با مجموعه‌ای از افراد. امکان مطالعه 
رها و وی گی‌هاع اضای افراه ماوت را که رنه 
شناسایی و کشف راه حل‌های کارآمدتر می‌شود. فراهم 
می‌سازد. در این مطالعه. الگوریتم ژنتیک رشته‌های متناسب 
با ارزش را برمی‌گزیند و آن دسته از رشته‌هایی را که 
تناسب کمتری با جمعیت مورد بررسی دارند حذف 
می‌کنند. هر کدام افراد جمیعت که تقریب‌هایی از جواب 
نهایی می‌باشند به‌صورت رشته‌هایی از حروف یاادغام 


متداول‌ترین حالت نمایش با ارقام صفر و یک است. حالت- 
های ورن استفاده از سه رقم اعداد حقیقی و اعداد صحیح 
هم مورد استفاده قرار می‌گيرند. برای مثال یک کروموزوم با 


دو متغییر ۵ و 0با ساختار شکل ۱ نمایش داده می شود. 
شکل ۱ - نمایش یک کروموزوم با ارقام صفر و یک 


مش ت اوه با تون سکس ا وا انار ان 
باقیمانده نمایش داده شده است. این می تواند به علت سطح 
دقت و يا محدوده متغییر تصمیم‌گیری باشد. مقادیر موجود 
روی کروموزوم‌ها به‌تنهایی معنی خاص ندارند بلکه باید از 
حالت کد شده خارج شوند تا به‌عنوان متغییرهای تصمیم- 
گیری دارای معنی و نتیجه باشند. باید توجه داشت که 
فرآیند جستجو بر روی اطلاعات کد شده انجام می‌گیرد مگر 
در صورتی که از ژن‌هایی با مقادیر حقیقی شود. بعد از اينکه 
کروموزوم‌ها از حالت کدگذاری شده خارج شدند می‌توان 
کارایی یا برازش هر فرد از جمیعت را محاسبه کرد. برازش 
انا ی اس گرا یشک افراه سای تلد مسا تست 


را نشان می‌دهد. در طبیعت برازش معادل توانایی فرد برای 


بقا می‌باشد. تابع هدف در تعیین برازش افراد نقش تعیین 
کننده دارد. در هنگام تکثیر به کمک اطلاعات اولیه‌ای که از 
تابع هدف به دست می‌آید. برازش هر فرد مشسخص می- 
گردد. از این مقادیر در فرآیند انتخاب استفاده می‌شود تا آنرا 
به‌سمت انتخاب افراد مناسب سوق دهد. هر چه برازش فرد 
نسبت به جمیعت بالاتر باشد احتمال بیشتری دارد که 
انتخاب شود. هر چه برازش نسبی آن کمتر باشد احتمال 
انتخاب آن برای تولید نسل بعدی کمتر می‌شود. وقتی که 
برازش تمام افراد جمیعت مشخص شد. هر کدام با احتمالی 
که متناسب با میزان برازش آنها است می‌توانند برای تولید 
نسل بعد انتخاب شوند. عمل تکثیر در الگوریتم ژنتیک برای 
رد و بدل اطلاعات ژنتیکی بین یک جفت يا تعداد بیشتری 


از افراد به‌کار می‌رود. ساده‌ترین نوع تکثیر تقاطع یک نقطه 


است دو رشته شکل را در نظر بگیرید اگر یک عدد صحیح 
از یک تا تعداد ارقام رشته منهای یک انتخاب شود و 
اطلاعات دو رشته را در دو طرف این دو نقطه عوض کنید 
دو رشته جدید به وجود می آیند که آنها را فرزند می‌ناميم 
به‌عنوان مثال اگر عدد شش را برای دو رشته شکل انتخاب 
کنیم نتیجه تقاطع یک نقطه‌ای به‌صورت شکل در می‌آید. 


شکل ۲- دو کرموزوم بعد از تقاطع (فرزندان) دو کروموزم 
قبل از تقاطع (والدین) 


این عملگر الزاما بر تمامی رشته‌های یک جمیعت اعمبل 
نمی‌شود بلکه برای اعمال آن بر یک جفت رشته یک 
اختمال نسیت: داده: هی شنود. بعند از این مرحعالته نا اعتب ال 
جدید عملگر جهش - روی رشته‌های تولید شده اعمال 
می‌گردد. در جهش, هر فرد به تنهایی با توجه به قوانیین 
احتمال می‌تواند تغییر کند. 

در نمایش دودویی رشته‌ها. جهش به معنای تغییر مقدار 
یکی از خانه‌های رشته از صفر به یک و يا از یک به صفر 
می‌باشد. به‌عنوان مثال جهش در هفتمین خانه اولین فرزند 
تولید شده در مثال قبل منجر به ایجاد رشته شکل ۲ می- 
گردد. 


شکل ۳- کروموزوم بعد از جهش 
بعد از مراحل تکثیر و جهش کروموزوم‌ها از حالت کد شده 
خارج می‌شوند و مقدار تابع هدف هر یک محاسبه می‌شود 
سپس به هرکدام برازشی اختصاص داده می‌شود. حال اگر 
لازم باشد دوباره مراحل انتخاب و تکثیر و غیره انجام می- 
کفون نطو زو فا نان ی رود که زاین فیط 
جمعیت جواب ها افزایش یابد. الگوریتم وقتی پایان می‌یابد 
که هدف خاصی برآورده شود. به‌عنوان مثال تعداد مشخصی 
تا اشیاه قرا مستای فا یفام ناگوان 


برسیم. 


روش ترکیبی سیستم‌های استنتاج فازی و الگوریتم ژنتیک 
شکل ۴ فلوچارت روش سیستم استنتاج فازی عصبی 
تطبیقی را نشان می‌دهد. پس از فراخوانی» مجموعه داده‌ها به 
دو قسمت آموزش و تست دسته‌بندی می‌شوند. که در این 
مطلعه داده‌ها به دو دسته ۵۰ و ۵۰ تقسیم می‌شوند و به- 
ترتیب جهت آموزش و تست مدل. مورد استفاده قرار می- 
گیرند. پس از دسته‌بندی‌هاءتولید سیستم استنتاج فازی انجام 
می‌شود. جهت انجام این کار دو روش پارتیشن بندی شبکه- 
ای" و خوش‌بندی تفریقی آوجود دارد. در این مطالعه با 
توجه به عملکرد بهتر پارتیشن‌بندی شبکه‌ای در مطالعات 
گذشتکگان (ابتهاج و بنک‌داری ۲۰۱۴ [۱۷]) از ایین روش 
استفاده می‌شود. در این مرحله همچنین تعداد عضویت و 
نوع توابع عضویت ورودی و خروجی باید مشخص شود. 
در این مطالعه از تابع عضویت گوسی استفاده شده است. 
تعداد توابع عضویت نیز که به‌صورت سعی و خطا تعیین 
شده است. برابر با ۳ در نظر گرفته شده است. پس از تولید 
سیستم استنتاج فازی الگوریتم آموزش شبکه باید تعیین 
شود. عموما جهت آموزش سیستم استتاج فازی عصبی 
تطبیقی از دو الگوریتم هیبرید و پس انتشار که الگوریتم‌هایی 
شناخته شده هستند. استفاده می‌شود. مطالعات اخیر در زمینه 
انتقال رسوب در مجاری فاضلاب‌رو (بتهاج و بنکداری 
۴ شان دهنده نشان دهنده عملکرد بهتر الگوریتم 
هیبرید نسبت به پس انتشار است. از این‌رو در این مطالعه 
عملکرد الگوریتم هیبرید در مقایسه با الگوریتم زنتیک به- 
عنوان یک الگوریتم تکاملی مشهور که توسط محققین برای 
مطالعات آینده پيشنهاد شده است. مورد ارزیابی قرا می‌گیرد. 
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پیش‌بینی با استفاده از داده‌های تست مورد ارزیابی قرار می - 
گیرد و اگر نتایج رضایت بخشی حاصل شده باشد. روند 
مدل‌سازی خاتمه می‌یابد وگرنه این روند از مرحله تولید 
سیستم استنتاج فازی, تا رسیدن به جواب مورد قبول دوباره 


تکرار می‌شود. 


مدل آزمایشگاهی 

در این مطالعه برای صحت سنجی مدل‌های عددی از اندازه- 
گیری‌های آزمایشگاهی توسط آب قانی (۱۹۹۳) اتاو 
همکاران (۱۹۹۹) [۱۹] و وونگ ویسس امجای و همکاران 
(۲۰۱۰) استفاده می‌شود. آنها مقادیر ۳۲ عدد فرود جریان, 
۷) غلظت حجمی رسوبات. 0/8 نسبت قطر متوسط ذرات 
رسوب به شعاع هیدرولیکی. 17/۸ نسبت قطر لوله (کانال) 
به سطح مقطع جریان. 1/0 نسبت شعاع هیدرولیکی به قطر 
لوله م۲ عدد ذره بی‌بعد (۷3ظ050805-1(/۷ 40۲ تسبت 
قطر متوسط ذارت به قطر لوله و ,۸ ضریب مقاومت کلی 
بار رسوب (ضریب مقاومت جریان رسوبی) را اندازه‌گیری 
نمودند. به‌بیان دیگر آنها تاثیر پارامترهای مذکور را روی 
تغییرات عدد فرود (۳) مورد بررسی قرار دادند. در جدول 
۱ مقادیر حداقل» حداکثر و متوسط مقادیر آزمایشگاهی آب 
قانی (۱۹۹۳ اتا و همکاران (۱۹۹۹) و وونگ ویسس 


امجای و همکاران (۲۰۱۰) مرتب شده تخت 


شکل ۴ - فلوچارت روش مورد استفاده در اين مطالعه 


جدول ۱- دامنه مقادیر آزمایشگاهی آب قانی (۱۹۹۳) اتا و همکاران (۱۹۹۹) و وونگ ویسس امجای و همکاران (۲۰۱۰) 


یط ی پارامتر 
1809۳701 9522۳401 4.678۲1 (/1) دبی 
۳/00 4.500۳2 1.00۳02 (ص) قطر لوله 
1.103۳-1 2.291۳1 2.000۳2 نی بان 
5.725۳2 1.138۳-1 1.200۳2 ( شعاع مزر سکن 
2.523۳2 8.136۳2 1.401۳-3 2 سطح مقطع 
6.675۳1 1.216۳4-00 2.370۳-1 (۳/9) سرعت 
00 9.300۳0 2.000۳1 (س) قطر متوسط ذرات 
3.017۳-3 6.00013 5.0704 تست کایال 
2.697۳2 5.320۳-2 1.290۳2 غتر فافع اوقت کی بان ارتتو نب 
1907۳-۷02 03-+1.280۳ 00-+1.000۳۷ غلظت بخجمی زسویاث 
امبووج6 7 حمبوووم2..... ممپوروو5.0 تا 
029]0.+4 1126۳401 00-+1.519۳ 


نحوه ترکیب پارامترهای ورودی 

همان‌گونه که بیان شد در این مطالعه با استفاده از مدل 
ترکیبی انفیس- الگوریتم ژنتیک مکانیزم انتقال رسوب مدل 
سازی می‌شود. برای رسیدت به این هدف مقدار عدد فرود 
جریان شبیه‌سازی می‌شود. برای حساسیت سنجی 
پارامترهای ورودی که شامل ۷ 4/۴ 92/۸ طل بیط 
0 و ,2 ۱۲۷ مدل‌های ترکیبی مختلف معرفی می‌شود. 
مدل‌های مذکور به هفت دسته که شامل دسته اول با یک 
پارامتر ورودی» دسته دوم با دو پارامتر ورودی, دسته سوم 
به سه پارامتر ورودی. دسته چهارم به ۴ پارامتر ورودی؛ 
دسته پنجم با ۵ پارامتر ورودی» دسته ششم با ۶ پارامتر 
ورودی و دسته هفتم با ۷ پارامتر ورودی مقادیر عدد فرود را 
مدل سازی می‌نمایند. ترکیب پارامترهای ورودی در شکل ۵ 


نشان داده شده است. لازم به ذکر است که در مطالعه حاضر 


عدد فرود 


۰ درصد اندازه‌گیری‌های آزمایشگاهی برای آموزش و ۵۰ 
درصد باقیمانده نیز برای تست مدل‌های عددی استفاده شده 


اشتت اه 


شکل ۵- نحوه ترکیب پارامترهای ورودی برای مدل‌های 
مختلف ۲۳۲۹-0۸( 


شاخص‌های آماری مورد استفاده برای بررسی دفت مدل- 
های عددی 

در این تحقیق به منظور ارزیابی دقت مدل‌های عددی از 
شاخحص‌های آماری ضریب تبیین (472 درصد میانگین مطلق 
خطا ۷۸۳ خطای جذر میانگین مربعات (5) و 
شاخص پراکندگی (50) که در قالب معادلات (۱۱) تا (۱۴) 


معرفی شده‌اند: 


2 


(۱۱) 


۹9 


(۱۳ 


(۱۳ 


در معادلات مذکور مقادیر ز(هسدوی) 61 
ز(فمه‌نممجه () رزق و ه 0( و 8 به‌ترتیب برابر 
مقادیر آزمایشگاهی, نتایج پیش‌بینی شده توسط مدل‌های 
عددی. میانگین مقادیر آزمایشگاهی و تعداد اندازه‌گیری- 


های آزمایشگاهی می‌باشند. 


بحث و نتایج 

مدل‌های برتر 

همان‌گونه که در بخش قبل عنوان شد در این مطالعه 
هفت مدل ۸۳۲۹-0۸ برای شبیه‌سازی عدد فرود 
جریان سه فازی تعریف شده است (مدل‌های شماره یک 
تا هفت). بر اساس نتایج مدل‌سازی. مدل شماره شش 
به‌عنوان مدل برتری که مقادیر عدد فرود را با یک 
پارامتر ورودی تخمین می‌زند معرفی می‌شود. این مدل 
تامرر 0 ایک سای شوب شته رای آسد سل 
در هر دو حالت آموزش و تست به‌ترتیب مساوی 
۴ و ۰/۷۳۱ تخمین زده شده است. همچنین مقادیر 
1۸/5۶ شاخحص پراکندگی (81) در وضعیت 
کی اوه ی قکب ای ان ۱7۶۵۷۱۹/۱۳ و 
۸ محاسبه گردیده است. در ادامه به بررسی مدل 
برتری که مقادیر تابع هدف را با دو پارامتر ورودی 


پیش‌بینی می‌کند مورد ارزیابی فرار می‌گیرد. مدل‌های 


۱ )1 ز([ 0660 1۳ 


)1 ۳ کر 1 / | 12 / 3 2 > ی 


11 


7 
1 ِ م00) من 1 1زا 0 


۳ 
"9 


1 11 
1۸۸/۲ 2-2 


۷ آ<ز 


نون 5 | 


1 1 
سس رسمار 7 ۱ - 7/55 
1۳/5 


(56۳۷۵۵ بن() 


شماره ۸ تا ۲۸ مقادیر عدد فرود را با دو متغیر ورودی 


شبیه‌سازی می‌کنند. با توجه به بررسی نتایج مدل‌هایی 
که مقادیر عدد فرود را با دو پارامتر ورودی تخمین می- 
زنند. مدل شماره ۸ از دقت بالاتری برخوردار است.. به- 
عنوان مثال در حالت آموزش برای این مدل مقادیر 
۲ و ۸۸5 به‌ترتیب مساوی ۵/۵۲۹ و ۰/۳۱۵ 
تخمین زده شده ابیتتا: این در حالی است که شاخحص 
پراکندگی در شرایط آموزش و تست مدل شماره ۸ به- 
ترئیب مساوی ۰,۱۷۶ و ۱/۰۸۷ بداست آمده شنت 
همچنین مقادیر ضریب تبیین در وضعیت تست این مدل 
برتری که مقادیر تابع هدف. که در اینجا عدد فرود 
جریان سه فازی است. را با سه متغییر ورودی پیش‌بینی 
های شماره ۲۹ تا ۶۳ با سه پارامتر ورودی تعریف شده- 
اند. در این بین» مدل شماره ۱ مدل برتر شناسایی تال 
مقادیر 7۳ برای مدل شماره ۱ در هر دو شرایط آموزش 
و نست به‌ترتیب مساوی ۷ و ۰/۹۶۲ نخمین زده 
شده است. با توجه به نتایج مدل مذکور. مقدار 
شاخص‌های آماری ۰1/۸۵۴ 2۸/5۳ و آک در شرایط 
نست به‌ترئیب مساوی ۷/۸۹ ۹ ۰*۴۵ و ۰/۱۰۲ بدست 


آمده است. همان گونه که در قسمت‌های قبلی عنوان 


شد. مدل‌های شماره ۶۳ تا ۹۸ مقادیر عدد فرود جریان 
سه فازی را به چهار چارامتر ورودی تخمین می‌زنند. با 
تجزیه و تحلیل نتایج مدل‌های مذکون مدل شماره ۶۵ 
مقادیر تابع هدف را با دقت بیشتری شبیه‌سازی کرد. به- 
عنوان مثال مقدار ضریب تبیین این مدل در شرایط 
آموزشق بت هدر تیب ستاو ۰/۹۴۳۸ و ۱/۹۲۷ 
محاسبه شد. همچنین مقادیر ]۸/5 و 97 برای مدل 
شماره ۶۵ در شرایط تست به‌ترتیب برابر ۱/۵۱۰ و 
۹ محاسبه شده است. علاوه براین مقدار 1/۸۳ 
حالت آموزش این مدل مساوی ۹/۹۴۵ تخمین زده شده 
است. در این قسمت به ارزیابی مدل‌های با پنج پارامتر 
ورودی پرداخته می‌شود. بر اساس نتایج شبیه‌سازی؛ 
مدل‌های شماره ۹٩‏ تا ۱۱۹ مقادیر تابع هدف را با پنج 
پارامتر ورودی تخمین می‌زنند. با بررسی کلیه این مدل- 
های مشخص شد که مدل شماره ۱۰۰ دارای دقت 
پیشتر و خطای کمتری است. بر اساس نتایج مدل شماره 
۰ مقدار شاخحص پراکندگی و ضریب تبیین در حالت 
آموزش به‌ترتیب مساوی ۰/۱۱۶ و ۰/۹۳۸ می‌باشد. 
همچنین در وضعیت تست این مدل مقادیر ۸۸۸ و 
به‌ترتیب برابر ۱۱/۹۸۰ و ۰/۵۷۹ است. همچنین 
هفت مدل متمایز ۸۳15-0۸ مقادیر عدد فرود را با 
شش متغیر ورودی مدل‌سازی می‌کنند. به بیان دیگر 
هي فا ۱ ای از شیارا 
ورودی هستند. با ارزیابی مدل‌های مذکور مشخص شد 
که مدل شماره ۱۲۰ دارای همبستگی بالاتری در مقایسه 
پااشای ال ها اس تب اق ال داد رون ال تاه 
های آموزش و تست این مدل به‌ترتیب مساوی ۰/۹۰۶ 
و ۰/۹۰۰ محاسبه شده است. همچنین مقادیر 1/۸۳ 
,۳ و 57 در شرایط تست مدل شماره ۱۲۰ به‌ترتیب 


مساوی ۸ + ,+ و ۰/۱۶۳ تخمین زده شده 


است. علاوه بر این مقدار شاخض پراکندگی در حالت 
آموزش این مدل مساوی ۰/۱۴۳ محاسبه شده است. در 
بخش بعدی مدلی که با کلیه پارامترهای ورودی مقادیر 
تابع هدف را شبیه‌سازی می‌کند مورد مطالعه قرار می- 
گیرد. مدل شماره ۱۲۷ تابعی از کلیه پارامترهای ورودی 
است. بر اساس نتایج مدل‌سازی توسط مدل شماره 
۷ نو وضنعیت: آمتوژشی مفتار فحریت. تسیز بان 
۱ است. همچنین برای اين مدل در حالت آموزش 
مقادیر ۸/۸۵۴ و 57 به‌ترتیب مساوی ۱۱/۳۳۸ و ۰/۱۶۶ 
محاسبه شده‌اند. در مقابل در شرایط تست مقادیر ۰1 
و ۸/58 به‌ترتیب مساوی ۰/۸۴۲ ۱۲/۷۱۲ و 
۰ بدشست: آمده استستا 

در این بخش از مطالعه. به ارزیابی توزیع خطای مدل- 
های برتر که توسط الگوریتم هیبریدی ۸۳15-0۸ 
شبیه‌سازی شدند پرداخته می‌شود. در شکل ۶ توزیع 
خطای مدل‌های مذکور در هر دو حالت آموزش و تست 
به‌تصویر کشیده شده است. به‌عنوان مثال در شرایط 
آموزش» برای مدل شماره ۶ تقریبا ۱۵ درصد نتایج 
مدل‌سازی دارای خطای کمتر از ۴ درصد به خود 
اختصاص داده‌اند. همچنین برای این مدل حدودا ۳۵ 
درصد عدد فرودهای شبیه‌سازی شده خطایی کمتر از 
۰ درصد دارند. علاوه بر این در شرایط آموزش مدل 
شماره ۸ تقریبا ۹٩‏ درصد نتایج این مدل خحطایی کمتر 
از ۲۰ درصد هستند. در حالت آموزش برای مدل شماره 
۱ نیز تقریبا ۶۲ درصد نتایج دارای خحطایی کمتر از ۸ 
درصد دارند. حدودا ۸۰ درصد نتایج مدل‌سازی شده 
توسط مدل شماره ۶۵ خطایی کمتر از ۱۶ درصد به خود 
اختصاص داده‌اند. برای مدل شماره ۱۰۰ نیز در شرایط 
آموزش ۳ درصد عدد فرودهای شبیه‌سازی شده 


خطایی کمتر از ۴ درصد دارند. همچنین برای مدل 


شماره ۱۲۰ در وضعیت آموزش حدودا ۷۰ درصد نتایج 
دارای حطایی کفتر از ۱۴ دراصد هستند: این دز خالی 
است که برای مدل شماره ۱۲۷ حدودا ۳۴ درصد عدد 
فرودهای شبیه‌سازی شده در وضعیت آموزش دارای 
درصد خطایی کمتر از ۶ درصد هستند. در شرایط تست 
نیز روند مذکور مشابه شرایط آموزش مدل‌های برتر 
است. به‌عنوان مثال. مدل شماره ۶ با خطایی کمتر از ۶ 
درصد حدودا ۲ درصد نتایج را تقریب زده است. 
همچنین در شرایط تست مدل شماره ۸ تقریبا ۹۶ درصد 
نتایج مدل‌سازی شده خطایی کمتر از ۱۸ درصد هستند. 
تقریبا ۷۷ درصد نتایج مدل شماره ۱ در وضعیت تست 
مقدار خطایی کمتر از ۱۲ درصد هستند. این در حالی 
است که حدودا ۵۱ درصد عدد فرودهای مدل‌سازی 
ات قبط مک تیار عم هار تا که نت۸ 
درصد هستند. برای مدل شماره ۱۰۰ نیز حدودا ۲۴ 
درصد نتایج شبیه‌سازی دارای خطایی کمتر از ۶ درصد 
هستند. علاوه براین حدودا ۷۵ درصد نتایج مدل شماره 
۰ خطایی کمتر از ۱۶ به خود اختصاص داده‌اند. 
همچنین در وضعیت تست مدل شماره ۰۱۲۷ حدودا ۶۹ 
درصد نتایج این مدل خطایی کمتر از ۱۴ دارند. 

بر اساس نتایج مدل‌سازی عدد فرود جریان سه‌فازی 
توسط مدل فرا ابتکاری ۸۳15-0۸ در بین ۱۲۷ مدل 
توسعه داده شده. مدل شماره ۸ به‌عنوان مدل برتر 
معرفی می‌شود. این مدل مقادیرعدد فرود را با دو پارامتر 
ورودی تخمین می‌زند. مقادیر عدد فرود شبیه‌سازی شده 
و مقایسه آن با مقادیر مشاهداتی نیز در شکل ۷ قابل 


مشاهده است. 


شکل ۶- نمودارهای توزیع خحطای مدل‌های برتر در 
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شکل ۷- نمودارهای پراکندگی در هر دو حالت آموزش 


(17210) و تست (1690) مدل هشت 


نتیجه گیری کلی 

اگر جریان مجرای فاضلاب نیمه‌پر باشد و جربانی که 
انتقال رسوب در آن به شکل بار بستر است وارد این 
جریان شود ولی سرعت جریان به‌حدی نباشد که مانع از 
ته تسش قنوافه مه ری مکی راهن تاه اس 
موضوع مقاومت بستر را به‌دلیل افزایش عمق بستر و 
کاهش سرعت. افزایش می‌دهد. ممکن است فرض شود 
که کاهش سرعت دلیلی بر کاهش ظرفیت انتقال رسوب 
جریان می‌شود که باعث ته‌نشینی بیشتر و احتمالا انسداد 
می‌شود. هر چند که در حقیقت شواهد آزمایشگاهی 
نشان می دهد که حضور بستر ته‌نشین شده در واقع به 
جریان اجازه می‌دهد که ظرفیت انتقال رسوب بیشتری 
به‌صورت بار بستر را داشته باشد. در این مطالعه با 
استفاده از یک روش فرا ابتکاری که با ترکیب شدن مدل 
انفیس و الگوریتم ژنتیک توسعه داده شد. جریان سه 
فازی درون مجاری فاضلاب‌رو مدل‌سازی گردید. برای 
مدلهای با یک پارامتر ورودی» مدلی که تابعی از 
نسبت قطر متوسط ذارت به قطر لوله (0/1) بود. عدد 
فرود جریان سه فازی را با دفت بالاتری تخمین زد. 
ال نوی مه تست شام هتای ور 
به‌ترتیب مساوی ۱/۰۵۷ و ۱۹/۱۳۴ محاسبه شد. 
همچنین برای مدل‌های با دو پارامتر 0۷ و 4/8 تخمین 
ژده برای اتن ملل در شالت آوزش برای: این سل 
مقادیر ۸/۸۴ و 5 به‌ترتیب مساوی ۵/۵۲۹ و 
۵ بود. برای مدل برتری که عدد فرود جریان را با 
سه پارامتر ورودی تخمین زد. مقادیر ۰۸۸۳ ۸5۲ 


و ]5 در حالت تست به ترئیب مساوی ۵ ۰/۸۴۰۳ و 


ورودی. مقادیر تابع هدف را تخمین زد. مقادیر 5 مقدار 2 در حالت‌های آموزش و تست این مدل به- 
و 27 برای مدل شماره ۵ در شرایط تست به‌ترتیب ترئیب مساوی ۶ و ۰۰۰ محاسبه شدند. بر اساس 
برتر با ۵ متغیر ورودی. در وضعیت تست این مدل پارامتر ورودی مشخصات جریان سه فازی را مدل‌سازی 
7/۸ بود. علاوه بر این مدل برتری که مقدار عدد ۴ داشت. 
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شکل ۱ - نمایش یک کروموزوم با ارقام صفر و یک 


۱ ۱ 
۱0 ۱۱ 


شکل ۲- دو کرموزوم بعد از تقاطع (فرزندان) دو کروموزم قبل از تقاطع (والدین) 


1 سس 1110010110110 
شکل ۳- کروموزوم بعد از جهش 


فراخوانی داده ها 


تولید داده های آموزش و تست - نوع تابع عضویت ورودی 


انتخاب موارد ورودی به شبکه: 


- تعداد تابع عضیت ورودی 


- نوع تابع عضویت خروجی 


تولید سیستم استنتاج فازی 


(پارتیشن بندی شبکه ای 
- هیبرید (پس انتشار خطا و حداقل مربعات) موش شک 
تب ش 
- الگوریتم ژنتیک ۱ 


ارزیابی مقدار پیش بینی شده رسم داده های پیش بینی برای: 
(:1) - داده های آموزش 


-داده های تست 


- انتخاب روش بهینه سازی 
- تلورانس خطا 
- تعداد تکرارها 


شکل ۴ - فلوچارت روش مورد استفاده در این مطالعه 
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شکل ۵- نحوه ترکیب پارامترهای ورودی برای مدل‌های مختلف ۸۳15-0۸ 
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درصد خطای تجمعی 


آستانه درصد خطای نسبی 
9۵ 9 < 
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حالت آموزش ۱0 
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مقادیر عدد فرود مشاهداتی 
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